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In quenched lattice QCD with standard Wilson fermions the quark prop- 
agator is computed very close to the chiral limit in the zero-temperature 
case. Starting from our experience with lattice QED we employ a modified 
statistical method in order to estimate reliably hadron masses. 

1 Introduction 

Chiral symmetry has remained a problematic topic in lattice gauge theories 
over the years. As is well-known, Wilson fermions explicitly violate the 
chiral symmetry, which can be hopefully restored in the limit k — > k c {(5) {(3 
and k denote the gauge coupling and hopping parameter, respectively). The 
QCD study at the experimentally established m^/rrip value is obstructed 
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near the chiral 
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by a huge amount of required computer resources and - particularly for the 
quenched approximation - by large fluctuations in the observables caused by 
near-to-zero eigenvalues of the fermion matrix ('exceptional configurations'). 
In |l| a variance reduction technique for estimating the pseudoscalar mass 
was proposed. In the case of QED the new estimator for m n was proved 
to work well close to n c , where the standard estimator fails. In this work we 
extend the study to the case of quenched lattice QCD. 

The Wilson action for lattice QCD is S = Sg + Sf, where Sq denotes the 
plaquette SU (3) gauge action and Sf the fermion part 

s F (u^^) = £ Ev^MWi , (i) 

/=l,2 x,y 

with the Wilson matrix M.(U) = 1 — kV(U), 

^xy = E [^V,x+fiPfi U xfl + 5y,x-jlP^Ui-^ (2) 

andP± = i±7M- 

The observables of interest are the 7r and p- meson correlators (masses). 
On every configuration {U} we determine the operators 

where t — — y^\. 

All simulations were performed on a 16 3 x 32 lattice at (5 = 6.0. The 
number of the gauge field configurations is 0(100). k- values were chosen 
between 0.1558 and 0.1570. The inversion of M. was done with BiCGstabI 
with restarting. 

2 Correlators and Masses 

The standard definition of the effective mass m e ff(r) is given by the equation 
cosh(m e// (r)(r + 1 - ^)) _ (r [M (r + 1)) 



cosh(m e// (r)(r-f )) (T M (t)) 



(4) 
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Figure 1: Scatter plots for the pion correlator. 



The plateau in the r-dependence of m e //(r) is supposed to define the 'true' 
mass m^p]. However, for k — > n c {(3) the values of T n (r) become strongly fluc- 
tuating. There is no reliable estimate for the averages (T n (r)), the approach 
to K c {0) seems extremely difficult |2|, [3j. 

In it was shown that a new estimator for the pseudoscalar mass can 
be defined which is identical to the standard estimator in the case of linear 
correlations 

y(x) = C • x, (5) 

where x = T n (r); y = T n (r + fc); > 0, and r ?r (r + fc), IV (r) are calculated 
on individual configurations. y(x) stands for the conditional average of y(x) 
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Figure 2: over r for two values of k. (o) denote the standard and (x) 
the improved estimator. 



at a fixed value of x. If (Q) holds with x ^ , then 



CVCt + AQ) _ (y) 

M _ <x> 



T^r + fc) ' 



(6) 



and the effective masses are determined by the r.h.s. of (0). The important 
observation |J is that the ratio r 7r (r + k)/T n (r) calculated on individual 
configurations does not suffer from near-to-zero eigenvalues. The average of 
this ratio is statistically well-behaved, in contrast to the ratio of averages. 

A typical situation is shown in Fig.|l| where r 7r (r+l) over T n (r) is depicted 
for r = 9 at several k's (note the different scales). We do not find any 
qualitative difference to the quenched compact QED 4 case. Therefore, we 
expect the improved estimator to be applicable also for quenched QCD. 
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Figure 3: m\ as a function of 1/ac. 



Fig.^| represents the effective mass m^-(r) for two different values of k. 
For the smaller k we obtain very good agreement between the two estimators, 
which with more statistics, should become identical. In the lower plot, near 
to k c , the standard estimator and its error bars cannot be trusted anymore. 
In contrast, the improved estimator gives a very nice plateau in miff at r > 5 
which still is statistically well-defined. 

The extrapolation m^ — > is shown in Fig.|^. The PCAC-like dependence 

is nicely fulfilled. In this figure we included also data from the QCDSF 
collaboration on 16 3 x 32 and 24 3 x 32 lattices Q for comparison. The data 
look well-consistent for coinciding lattice sizes. 
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Figure 4: Scatter plot for p-meson correlator. 



The p-meson mass m p should not disappear in the limit k — > k c . There- 
fore, the p-meson correlators V p {r) are expected to be less sensitive with 
respect to near-to-zero eigenvalues. The standard estimator for m p works 
better, in contrast to the pion case. However, it is interesting to note that 
scatter plots for the p-meson correlators with different r show again the 
linear correlations between different r-slices (see Fig. 
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3 Summary 



We observed that due to strong linear correlations between the values of the 
pseudoscalar correlators at different r the ratios r 7r (r + l)/r 7r (r) do not suffer 
from near-to-zero (exceptional) eigenmodes of the fermionic matrix. This 
means that for every given configuration all 'divergent' contributions to the 
correlators are factorized. 

Making use of this observation we proposed another estimator of the pseu- 
doscalar mass, which is well-defined near K c (f3) in contrast to the standard 
estimator. 

For K-values sufficiently below k c (0), i.e., where the standard estimator 
can be reliably defined, both estimators are shown to be in a very good 
agreement. By approaching k c (/3) we observe a PCAC-like dependence as in 
eq. (0). For values of k very close to k c (/3) the standard estimator fails to 
work, while the improved one still fits the same straight line having a very 
small statistical error. 

The new estimator permits to approach the chiral limit much closer than 
the standard one would allow. Especially, systematic errors induced by the 
finite volume and/or by the quenched approximation can be investigated 
in more detail within the 'critical' region. In particular, the 'critical' value 
k c (P) can be determined with higher accuracy. 

An alternative way to cure the problem of 'exceptional configurations' 
was proposed in Ref. ||. 
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